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Abstract 

After a short review of various ways to calculate the Maslov index appearing in semiclassical 
Gutzwiller type trace formulae, we discuss a coordinate-independent and canonically invariant for- 
mulation recently proposed by A Sugita (2000, 2001). We give explicit formulae for its ingredients 
] and test them numerically for periodic orbits in several Hamiltonian systems with mixed dynamics. 

We demonstrate how the Maslov indices and their ingredients can be useful in the classification of 
periodic orbits in complicated bifurcation scenarios, for instance in a novel sequence of seven orbits 
born out of a tangent bifurcation in the Henon-Heiles system. 
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1 Introduction 

The idea of Gutzwiller P to express the quantum-mechanical density of states of a Hamil- 
Q ! tonian system in terms of its classical periodic orbits and their properties has brought an 
| enormous progress to the field of "quantum chaos" [2J. The trace formula which he derived 
in 1971 pQ is valid for systems in which all periodic orbits are isolated in phase space. Later 
versions of trace formulae were formulated for integrable systems [31 E] and mixed systems 
with continuous symmetries [31 El E] , arid very recently also for systems with explicit spin 
degrees of freedom E] . 

An important ingredient in all these trace formulae is, besides the action and stability of a 
periodic orbit, the so-called Maslov index. It is an invariant property of a periodic orbit which 
can change only when the orbit undergoes a bifurcation or when a continuous symmetry is 
broken or restored under the variation of a system parameter (eg energy, deformation or 
an external field). The calculation of the Maslov index is not always straightforward, in 
particular for systems with many degrees of freedom or systems which are not of the "kinetic 
plus potential energy" type. In the standard methods used in the literature [01 EH ITT] , 
the determination of the Maslov index of a stable orbit necessitates the explicit use of an 
"intrinsic" coordinate system that follows the orbit (as introduced by Gutzwiller [JJ), which 
can be numerically quite cumbersome. Easy-to-use calculational recipes using the method of 
^ ! JJJ have been given in the appendix D of |12j . 

Recently, Sugita J31 has given a formula for the Maslov index which only contains canon- 
ically invariant ingredients. However, no practical recipes were given for the explicit calcu- 
lation of the winding number m [see equation (j7j) in section 12. 2j . Inspired by another recent 
work on semiclassical trace formulae and Maslov indices ^1] , we have developed simple calcu- 
lational recipes for the winding number and the stability angle used in Sugita's formula [T3j . 
The purpose of the present paper is to implement these recipes for various Hamiltonian sys- 
tems with mixed dynamics and to test their results towards those obtained with the standard 
methods [HEBE]- As an outcome, we shall give some empirical rules for the behaviour of 
some of the ingredients of the Maslov indices at bifurcations and demonstrate that they can 
be useful for the classification of periodic orbits in connection with complicated bifurcation 
scenarios. 
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2 General formalism 



2.1 Trace formulae and role of Maslov indices 

For isolated orbits, Gutzwiller's trace formula pQ for the oscillating part of the density of 
states g(E) reads 



1 _ „ ,„ « 1 



^ppo ,=i v /|det(M r ; d -/)| 



8g{E) = -L J2 T ppo {E) £ , / cos \ S ppo (E) - 1 



(1) 



(2) 



where T ppo and S ppo are the periods and actions of the primitive periodic orbits (ppo) and 
r is an index which counts the repetitions of the primitive orbits. M re d is their reduced 
monodromy matrix - often called stability matrix - and a r is the Maslov index of the repeated 
orbit (which for stable orbits is not a multiple of that of the primitive orbit). For isolated 
orbits the Maslov index is always an integer, irrespectively of the dimensionality of the system. 

For systems with continuous symmetries, periodic orbits are no longer isolated in phase 
space but come in degenerate families. The semiclassical trace formulae can still be written 
in the general form 

Sg(E) = £ A po (E) cos [~ S po (E) - - a, 
po Lft a 

where "po" now refers to all repeated families and isolated orbits. The amplitudes A po (E) 
depend on the degeneracies of the orbit families, besides their stabilities, and have been given 
in the literature jHUHUE]. The Maslov index o~ po of an orbit family can be a half integer. Its 
determination, especially for families with high degeneracies such as they occur in isotropic 
harmonic oscillators with SU(N) symmetry, is by no means trivial. Although some hints 
may be found in (HI El there exists, to our knowledge, no simple recipe for their 

calculation. We shall in the remainder of the present paper restrict ourselves to systems with 
isolated orbits. 

It has been shown ^TJ ^] that the Maslov index of an isolated orbit is a canonical 
invariant, independent of the orbit's intrinsic coordinate system and of the starting point 
for the calculation of the monodromy matrix. However, the published practical ways for 
its calculation [TJ UHl IH] do make explicit use of the orbit's intrinsic coordinate system. 
Furthermore, using the form 

Opo = /ipo + fpo (3) 

given in [TJ] - whereby the part /x po is the Morse index appearing in the semiclassical Green's 
function and the part v po comes from the stationary-phase integrations transverse to the orbit 
- it has been explicitly demonstrated for an unstable orbit that \x and v alone need not be 
invariant and may, in fact, depend explicitly on its starting point. 

A remark on the name "Maslov index" might be in place here. Originally, the Maslov 
index was introduced in the framework of the WKB and later the EBK quantization of 
integrable systems [T71 HH1 EH] . The quantity a po appearing in the above trace formulae is a 
different index, although for integrable systems it can be easily related [U IE] to the Maslov 
indices of the EBK quantization. In the standard literature on the periodic orbit theory, the 
name of a po has, however, established itself as "Maslov index" and we want here to adhere 
to this convention. 

The trace formulae (0) and (0) provide an approximative description of the quantum- 
mechanical density of states in terms of classical periodic orbits and their properties. It is 
obvious that the Maslov index a po plays a vital role in establishing the correct quantum phase 
interferences and therefore must be determined correctly. 
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2.2 Invariant calculation of the Maslov index 



All information about the Maslov index of an isolated orbit is contained in its matrizant 
M(t) describing the time propagation of a small perturbation Sq, Sp around the orbit in 
phase space: 

(2$) -"<«(£?>)• 

where q(t) and p(t) are ^-dimensional coordinate and momentum variables. M(t) is ob- 
tained by solving the linearized equations of motion of a classical system characterized by its 
Hamiltonian H(q,p), which leads to the differential equations 

^M(t) = JH"(t)M(t), M(0) = I 2N , (5) 



where 



d 2 H d 2 H 
dqdq dqdp 



J =[-i„ 5 < s "c>= A* ■ < 6 > 

\ dpdq dpdp I 

J is the unit symplectic matrix in the 2iV-dimensional phase space, and /27V, In are the 2N- 
and iV-dimensional unit matrices, respectively. At the time of the orbit's period T, the matrix 
M(T) is called the monodromy matrix. One can always transform M(T) into block form. 
One parabolic (2x2) block contains the trivial unit eigenvalues related to energy conservation; 
the remaining (2N — 2)-dimensional part of M(T) is called the reduced monodromy matrix 
M re d or stability matrix. M re( j falls into (inverse) parabolic, elliptic, or (inverse) hyperbolic 
(2x2) blocks or, for N > 2 dimensions, loxodromic (4x4) blocks, depending on the stability 
of the respective orbits. 

The amplitude in the Gutzwiller trace formula (P) diverges when det(M r r ed — /) becomes 
zero. For isolated orbits this happens at bifurcations, where an orbit changes from stable to 
hyperbolically unstable or vice versa, or when a continuous symmetry is restored under the 
variation of a system parameter such as energy, deformation or some external field. At such 
singular points M r r ed contains an extra parabolic block. The trace formula (0) then cannot 
be used and one must resort to uniform approximations [201 HU 1221 I2H] which we, however, 
shall not be concerned with in the present paper. We just mention here that bifurcation and 
symmetry-restoring points are the only places where a Maslov index is allowed to change. 
The corresponding rules for bifurcations can be found in |2*T] : an example for the changes of 
Maslov indices under symmetry breaking (or restoring) will be mentioned in section 12.31 

Referring to the block form of the stability matrix, Sugita jT3| has recently proposed the 
following general formula for the Maslov index a r in (0) of an isolated periodic orbit: 



"ell / , 

<x r = Wl+2 /Vi 



i=l 



2vr 



+ rriih + 2mr . (7) 



Here [x] denotes the integer part of x (ie the largest integer number < x), n e n is the number 
of elliptic blocks and Xi the stability angle of the ith elliptic block, is the number of 
inverse-hyperbolic blocks, and m is a winding number. Sugita has also discussed in ^3] the 
relation of the winding number m to the homotopy theory. However, no explicit prescription 
for the computation of m has been given. 

This question has been partially answered in a recent review article on periodic orbit the- 
ory [Tl] . The winding number m has been identified as the Gel'fand-Lidski winding number 
and it has been explained how to calculate it in principle. Muratore-Ginanneschi 



3 



has also discussed a relation of a r given by (JZJ) to an index which is known in the mathe- 
matical literature after the names of Conley and Zehnder [22], and has provided extensive 
references on the latter subject. However, in ^3] the way of extracting a unique value of the 
stability angle Xi from the the eigenvalues e ±%Xi of the stability matrix has only been hinted 
at, and a practical algorithm still remained to be specified. We also want to remark that a 
representation similar to (J7J) appears in a mathematical paper where a classification of 
the admissible normal forms of the elements of Sp(2N) is given. 

The goal of the present section is to specify the definitions of Xi an d rn, in order to make 
the definition of a r in (J7|) unique and useful for practical applications. In particular, we are 
going to establish how the choice of Xi made in ^Hj corresponds to the prescription for the 
calculation of m given in [T4*] . 

Following [Hj, we split M(t) into a product of a periodic and an average part (also called 
the Floquet decomposition): 

M(t) = M per (t)M av (t) (8) 

with 

M w (t) = exp(tK) , (9) 

where K is a constant matrix. By definition, the periodic part of the matrizant in (jHJ) satisfies 
the condition M per (t) = M per (t + T). In particular, M per (0) = M per (T) = I 2 n- We therefore 
can specify the constant matrix K by equating 

M av (T)=exp(TK) = M(T). (10) 

Then, we get 

ff = ~ln[M(T)]. (11) 

To take the logarithm on the rhs of JHJ, we diagonalize M(T), calculate the logarithms of 
the eigenvalues of M(T), and then return to the initial basis. However, the relation fTTj) 
remains symbolic until we adopt a certain phase convention for determining the eigenvalues 
ofK. 

In the standard definition of the function ln(,2) with z = \z\e l< ^, the phase range <fi G (— ir, it] 
is chosen, corresponding to the branch cut line being taken along the negative real axis. Let 
us consider the eigenvalue problem 

M(T)if = e ± ^it, (12) 

with & = [£+]* and 

Xi = -iln[e +i »]e(0,7r). (13) 

The case Xi = 71 wm be discussed separately below. 
Let us now introduce the symplectic product 

Sl = +[Re(C + )] T JIm(C + ) = -[Re(C)] T ^Im(C) , (14) 

known as the Krein invariant PHI EH] • With this, we can adopt the following convention for 
the eigenvalues ±ijr of K: 

K£f = ±^if, (15) 



so that 



Xi = Xi for Si > 

Xi = 2n-Xi for ^ < . (16) 
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These relations fully determine the constant matrix K and specify uniquely the stability 
angle to be used in the formula (JJJ). In this phase convention, Xi takes values in the range 
(0,2vr). 

The case of an inverse parabolic block with e ±%Xl = — 1 is degenerate and requires special 
consideration. It occurs when trM re d = —2 and corresponds to the stability changing between 
elliptic and inverse hyperbolic. In this case we choose the value Xi = ^ by continuity reasons. 
The inverse parabolic block should be taken into accout in the formula (|7j) as a special case 
of either an inverse hyperbolic or an elliptic block, but not twice - in order to avoid double 
counting. 

The winding number m is an invariant characteristic of M peT (t) = M(t)M~ v l (t). To 
determine it, it is convenient to employ the so-called polar decomposition of the symplectic 
matrix M per into a product of an orthogonal matrix R peT and a positive-definite symmetric 
matrix W per : 

M per = R peT W pei . (17) 
In turn, the orthogonal matrix R pev admits the representation 



Rper 



Y V 

^-per 1 per 

^per A per 



Therefore, the winding number m can be defined as 

m = <p{T) - V (0) , (19) 

where 

<p(t) = — Argdet [X pev (t) + iY peI (t)} . (20) 

Since X per (t) and Y per (t) are periodic, m is a (positive or negative) integer number. 

The winding number has been vastly discussed in the literature, both mathematical 
and physical (see ^3] for extensive references). In particular, we would like to quote here 
that it has been introduced in [23] for a topological characterization of the structural stability 
of linear Hamiltonian flow. 

The extraction of R per (t) from M per (t) provides a nice representation of the evolution of 
<p(t), because det[X per (t) + iY per (t)} runs around the unit circle. However, the polar decom- 
position ()17j) is not essential for the calculation of the winding number m, even though the 
latter is encoded in R per (t). The same result as in (f27I|) can be also obtained from if)(T) —i/;(0), 
where 

i){t) = —Argdet [A pcv (t) + iB per {t)} , (21) 
and the matrices A per (t) and B peT (t) are the blocks of 

- ( few tl) ) ' (22) 

For a proof and further discussion of this point, see the Appendix A of [27] . 
Let us now consider the following canonical transformation 

M(t) = S , (t)M av (t)S'(0)- 1 . (23) 

with S(t) = M per (t). Since 5(0) = I 2 n, this expression is equivalent to (jHJ)- The relation 
between the Maslov indices a of M{t) and cr av of M av (t) for r = 1 is given by [IHj 

a = a av + 2m . (24) 
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But the winding number in cr av equals zero, since M av (t) belongs to the same homotopy class 
as the identity matrix, ie it can be continuously shrunk to the latter. Therefore, in order 
to determine a av we just need to find the number of elliptic and inverse hyperbolic blocks, 
which can be read off the block form of M av (T) = M(T). We remark that neither <7 av nor 
m depend on the choice of the starting point on the periodic orbit as it must be for 
canonically invariant quantities. 

We note in passing that the Maslov index a r can be identified with the winding num- 
ber obtained from a polar decomposition of the whole matrizant M(t) |TB]. However, that 
approach also requires some further specifications for stable orbits. 

This completes our specification of Sugita's approach. Before discussing another choice 
of the phase convention which allows us to make contact with the earlier approach of [TT] . 
we will illustrate our method with an analytical example. 



2.3 Analytical example: irrational harmonic oscillators 

We consider here a simple integrable system with isolated orbits, for which all the above 
quantitites can be evaluated analytically. This is the two-dimensional anisotropic harmonic 
oscillator 

H= 1 -(pI+P 2 v ) + 1 + uV) = £ (Pl + Ql) + ^ ( R 2 + Ql) , (25) 



where Q x = xJuT x , Q y = yJuJZ and P x 



Vvly/Uy- We assume that the 



Px/y/^x~, Py 

frequencies u x and u y are incommensurate, so that the only periodic orbits are librations 
along the x and y axes; they are isolated and stable. For the orbit along the x axis, the 
period is T x = 2tt/uj x , and the monodromy matrix and its periodic and average parts are, 
respectively, given by (cf [T2"l 



M x (t) 



M X;Pei {t) 



I 


cos(a; a; t) 










COS(L0yt) 




— sin(a; 3 .t) 


( 


V 





— sm(ujyt) 




cos(u x t) 










cos(Acjj / 1) 




— sin(a; a ;t) 





V 





— sin(Au;j / 1) 



sin^i) \ 

sm.{u)yt) 
cos{u x t) 

COS{<jJyt) ) 

sm{uj x t) 



where 



/ 1 



Vo 



— LOy 





cos{u x t) 


\ 

cos(a)yt) sm(uj y t) 

1 
- sin(a; y t) cos(u) y t) J 



sin(A ujyt) 


cos(Aa; ?; t) J 



AWy UJy UJy Ul X 



0J n 



In order to calculate the winding number m x , we consider 

Argdet [X x>peI (t) + iY XiPei (t)] = (uj x + Au y ) t . 

Then we obtain easily 

UJ X + A^y 

m x = = 1 + 



(26) 



(27) 



(28) 

(29) 

(30) 
(31) 
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Next we calculate Xx/27r: 



X.X 

2^ 



UJy 



UJy_ 



Finally, we obtain the Maslov index for the rth repetition to be 



1 + 2 



2tt 



+ 2rm 3 



1 + 2 





'UJy' 






lUJy- 




UJy 




+ 2r (l + 


J = l + 2r + 2 


- u x 












- Vx- 



(32) 



• (33) 



This result agrees with that obtained in [2B] using the method of [TT] . 

A note about the isotropic harmonic oscillator with uj x = uo y = uj may be of some interest 
here. The periodic orbits in this system, due to its SU(2) symmetry, are not isolated but 
form families of two-fold degenerate orbits. (The same is true for arbitrary rational axis ratios 
u x : Uy = n : p with integer n,p for which the orbits are Lissajous figures.) As mentioned 
in the introduction, we lack a general prescription for the calculation of the Maslov index 
of these families. However, the semiclassical trace formulae of isotropic harmonic oscillators 
can be obtained by other means and are found to be quantum- mechanically exact j2H]- In 
two dimensions, the trace formula reads 



9(E) 



E 



(huo) 



1 + 2 ^2 cos ! ' 

r=l 




(34) 



which supports a Maslov index ay = (mod 4). Indeed, with our above results we find 
m = 2 from ()31|) and \ = from (|32|). leading to = 4r which is equivalent to (mod 4). 
The reason for our identifying this Maslov index for the isotropic harmonic oscillator here 
is that in the context of perturbation theory, trace formulae for slightly perturbed harmonic 
oscillators have been developed |2~H] EHj in which the Maslov indices of the perturbed isolated 
orbits are obtained analytically, once the value for the unperturbed families is known. For 
instance, in the Henon-Heiles system discussed later in section E21 the unperturbed families 
break up into three isolated orbits A, B and C as soon as the nonlinearity is turned on [e > 
in ([44)1 ]. The changes in their Maslov indices with respect to af^ were found analytically 
to be Acta = +1, Acre = and Aa c — ~ 1- Indeed, the numerical methods for the isolated 
orbits yield oa = 5, Ob = 4 and oq — 3, both using the method of [TT] (see [HI]) and with 
our present method (see table EJ below). 

A straightforward generalization for the ^-dimensional harmonic oscillator with irrational 
frequency ratios uji/ujj = 1, 2, . . . , N) yields the Maslov index for the orbit along the j 
axis 



*£> = {N-l)+2r + 2Y, 



N 



i=l 



UJa 



(35) 



2.4 Alternative prescription and relation to earlier approaches 

In this section we introduce another prescription for calculating the quantitities \% an d m 
appearing in [IB]. It is based on the alternative Floquet decomposition 

(36) 



M(t) = Mper(t) M av (t) = M peI (t) exp(tK) 
specified by a constant matrix K such that 



(37) 
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This actually represents another convention for the choice of the stability angle. The formula 
for the Maslov index is then modified to 



"ell 



i=l 



£ 1+2 



sign (Si 



rXi 
2tt 



+ TUfa + 2 mr . 



(3f 



If Si > for all i, we have K = K and Xi 



Xi 



as well as m 



m. 



Then, there is no difference 



between and (|ZJ). If Sj < for some z, we can make the transformation 



rXi 
' 2tt 



-2r + 2 



2tt 



-2r + 2 



2tt 



(39) 



which is the winding 



Correspondingly, the winding number m of M per (t) changes to m, 
number of M per (t), such that 

2mr = 2mr + 2r . (40) 

Summing up (|3~9j) and (}40j) . we see that the o> in both (|3*8*jl and ((ZJ) coincide. Thus, the 
equivalence of both representations is established. We also note that, in general, the difference 
(m — m) equals to the winding number of e l ^ K ~ K ) ( which is the number of elliptic blocks of 
K (or K) with negative values of Sj. 

The sign of may change from positive to negative (or vice versa) away from bifurcation 
or symmetry restoring points. As a consequence, m changes its value by +1 or —1, but such 
as to conserve the total Maslov index. In two-dimensional systems, we have found this to 
happen when the stability discriminant 2 — det(M rec j — h) = trM rc d crosses or touches the 
line trM rec j = —2. The prescription for K based on (|3*Tj) is not relevant from the point of 
view of a canonically invariant formulation, but such a representation often appears to be 
more convenient in numerical computations. It reveals itself useful, in fact, to establish some 
relations to the approach of Creagh, Robbins and Littlejohn fl] (see also jTHI and Appendix 
D in J2])- A s mentioned already above, these authors have written the Maslov index as 
a sum (jSJ) of two contributions, which for stable orbits must be calculated separately. For 
stable orbits, \i and v are invariants in the sense that they do not depend on the starting 
point of the orbit, but they may change their values away from bifurcations or symmetry 
restoring points, exactly as it happens for m. We note that for a two-dimensional system, v 
is given by the upper right element b of the stability matrix 



M, 



red 




and can be calculated as 



dr ± {T) 



(41) 



(42) 



where p±(t) and r±(t) are the momentum and coordinate, respectively, transverse to the orbit. 
This relation actually demonstrates the necessity of knowing the orbit's intrinsic coordinate 
system for the calculation of b and hence of v. 

We conclude this section by a number of rules for the Maslov indices for two-dimensional 
Hamiltonian systems that are useful for the classification of periodic orbits, in particular in 
connection with complicated bifurcation scenarios such as will be discussed in our applications 
in sectional The rule 1 is rigorous and follows directly from the formulae (J7J) and (|38p. whereas 
the rules 2-4 are empirical, being based on numerical experience. They will be illustrated 
in the examples given in the next section. 
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1. In two-dimensional systems, the Maslov index oy is always even for hyperbolically 
unstable orbits and odd for stable and inverse-hyperbolically unstable orbits. 

2. For stable orbits, sign(s) = sign(5). This means that cr av is always negative when v = 
and positive when v = 1. 

3. The values of sign(s) of two stable orbits involved in the same bifurcation are identical. 

4. All orbits involved in a bifurcation locally have the same winding numbers m. 

These rules are consistent with the fact that Maslov indices of isolated orbits can only 
change at bifurcations and in symmetry-restoring limits. The changes at bifurcations are 
given in the papers of Sieber and Schomerus [21] and were found to be correctly reproduced 
by the present method in all cases. 

3 Numerical applications 

In this section, we shall apply our method to some systems with mixed classical dynamics. 
We first discuss two textbook systems, the homogeneous two-dimensional quartic oscillator 
and the famous Henon-Heiles system. They both have Hamiltonians of the form H = p 2 /2 + 
V(q) and have been investigated numerous times in the framework of periodic orbit theory 
using the previous methods [TU] and [TT] for the calculation of the Maslov indices. In the 
next two examples we shall study systems with spin degrees of freedom, for which a simple 
separation into kinetic and potential energy is not possible and the previous methods are not 
straightforwardly applied. 

3.1 The quartic oscillator 

As a typical system which exhibits the transition from integrable regular to almost completely 
chaotic motion, we study the quartic oscillator Hamiltonian 

H = \{vl + v 2 u ) + \{x' + y 4 ) + ^ xY- (43) 

It is homogeneous in coordinates and momenta, so that the energy can be scaled away. The 
chaoticity parameter is a. For a = 0, 1 and 3, the system is integrable and in the limits 
a —>■ — 1 and a — > oo it becomes nearly chaotic [32] • The stability of the linear orbit running 
along either of the axes, which we here denote by A, is known analytically |33; . The trace of its 
stability matrix is given by trM red = 4 cos(tc\/1 + 8a/ 2) + 2. Isochronous bifurcations of the 
primitive A orbit, which are of pitchfork type, occur when trM re( j takes the value +2 which 
occurs at the values a n = |n(n + l) with n = 0,3,4,5,... Period-doubling bifurcations 
of island-chain type occur at the values a p = 2p (p + 1) + 3/8 with p = 0, 1, 2, . . . These 
bifurcations of the A orbit and the analytical properties of the period-one and period-two 
orbits created at the bifurcations have been discussed in [31] and [3E]. As an example for 
our present evaluation of the winding number m, we shown in the left part of figure Q the 
function ip(t) defined in (j2"U|) . obtained here for the stable orbit L7 born at the bifurcation of 
the A orbit at a = 15. It yields m = 3. Since s > in this case, ip(t) is identical to (p(t). 
The Maslov index becomes <j\ = 7, in agreement with the value obtained in [31] using the 
formalism of jTTj. In the right part of figure ^ we show the shape of the L 7 orbit in the (x, y) 
plane. 
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t/T x 



Figure 1: Properties of the the L7 orbit in the quartic oscillator Hamiltonian Q43[) at a = 16.8. Left 
panel: phase function ip(t) (|20|) giving m = 3; right panel: shape of the L7 orbit in the (x, y) plane. 

In [HE] , the scaling behaviour of the fix points corresponding to period- four orbits created 
at bifurcations of the fourth repetition of the A orbit (denoted here by A 4 ) have been dis- 
cussed. To illustrate this scenario, we show in figure El the stability discriminant trM ret j of 
the period- four orbits involved in an island-chain bifurcation of the A 4 orbit, which occurs at 
a = 5.4305556, and in the succeeding pitchfork bifurcations of the P21 and P 21 orbits. The 
subscripts of the orbit names indicate their Maslov indices as obtained using the formulae of 
[TT] . (Not shown are the orbits created at the bifurcations of A 4 occurring at a = 4.375 and 
at a = 6.) 




a 

Figure 2: Stability discriminant trM re( j of the period-four orbits in the quartic oscillator (|43j) 
involved in an island-chain bifurcation of the A 4 ^ orbit occurring at a = 5.4305556 and in the 
succeeding pitchfork bifurcations of the P21 and P 21 orbits. 
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Figure 3: Poincare surface of section (y,p y ), shown in the left panel, taken at x = in the quartic 
oscillator (|43[) at a=6.0. The fix point in the central stability island corresponds to the A orbit, and 
the chain of eight unstable and stable fix points correspond to the two pairs of period-four orbits 
F22 and P21 with shapes shown in the upper and lower right panels, respectively. 

In figure El we show on the left side a Poincare surface of section (y,p y ), taken for x = 
at a = 6.0. The fix point in the central stability island corresponds to the A orbit and its 
repetitions. The KAM chain of eight unstable and stable fix points, which form the boundary 
of the stability island towards the chaotic sea, correspond to the two pairs of period-four orbits 
F 22 and P21 whose (x, y) shapes are shown in the upper and lower right panels of figure El 
respectively. These are the fix points whose scaling with a was studied in 

Tabled shows the stabilities, Maslov indices and their ingredients of the above period-four 
orbits in the quartic oscillator. We give the intervals (ct m i n , a max ) in which the orbits are 
stable (ell) with a fixed sign of s, hyperbolically (hyp) unstable, or inverse-hyperbolically 
(i-hyp) unstable. Values of (cc TO j n and a max ) marked by an asterisk (*) denote bifurcation 
points for the orbits listed in the corresponding rows. Note that the values of \x and v are 
not unique for unstable orbits; they may depend on the starting point along the orbit chosen 
for their calculation, but such that o = fi + v is invariant ^T] . For stable orbits, they are 
constant in each of the given regions and related to 2m and <r av as discussed in section 12.41 
Our results show the consistent agreement between the definitions of the Maslov index by 
Creagh et al [TT] and by Sugita [T3|. They also illustrate the empirical rules given at the end 
of section l2~4l 
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Table 1: Stabilities, Maslov indices and related properties of the orbits in the quartic 
oscillator whose stability discriminants tr M re d are shown in figure|2l 'ell', 'hyp' and 'i-hyp' 
denote elliptic (stable), hyperbolic and inverse-hyperbolic (unstable) orbits, respectively. 
Values of a m i n or a max marked by an asterisk (*) denote bifurcation points for the orbits 
given in the corresponding row. For unstable orbits, the decomposition of a = fi + v is 
not unique. 



3.2 The Henon-Heiles system 

Another famous system with mixed classical dynamics is given by the Henon-Heiles Hamil- 
tonian 

h = 2 (pl + p 2 y ) + 5 + y 2 ) + e ( x2 y ~\y 3 ^ ( 44 ) 

where e regulates the chaoticity of the system. The potential in f)44j) has three saddle points at 
the energy E* = 1/6 e 2 , over which a particle can escape if E > E*. The classical dynamics 
depend only on the scaled energy e = E/E* = 6e 2 E; in this variable the saddles are at 
e — 1. Along the symmetry lines passing through the saddles, one of them being the y axis, 
there are librating orbits (denoted here again by A) whose stability oscillates infinitely many 
times as the energy approaches the critical value e = 1, giving rise to an infinite cascade of 
isochronous pitchfork bifurcations. The scaled bifurcation energies e n (n = 1,2, ... , 00) form 
a sequence that cumulates at eoo = 1 in a Feigenbaum-like fashion; the orbits born at the 
bifurcations exhibit self-similarity with analytically known scaling constants [38 . 

At each successive bifurcation e n , the A orbit increases its Maslov index by one unit. The 
orbits born at the bifurcations are alternatingly stable rotations R CT and unstable librations 
L CT ; they can be uniquely classified by their increasing Maslov indices: R5, L 6 , R7, L 8 , etc, 
according to the rules given at the end of section 12.41 fcf also jHHEHl)- Besides the A orbit, the 
system possesses a curved librating orbit B which is unstable at all energies, and a rotating 
orbit C which is stable up to e = 0.8919 where it turns inverse hyperbolically unstable. Its 
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Table 2: Shortest period-one and period- two orbits in the Henon-Heiles system, their 
stabilities, Maslov indices and related properties. Notation as in tabled e = E/E* is the 
scaled energy; its values denoted by asterisks (*) are bifurcation energies e n . 

second repetition bifurcates at this energy, giving birth to an orbit D that stays stable up to 
e = 1.2375 where it becomes inverse hyperbolically unstable. We have calculated the Maslov 
index of all these orbits using the formulae given above and verified that they agree with the 
values obtained in [23 EU EH] using the method of [TT] and in [31] using the method of [TU] ■ 
The results are given in tabled again in energy intervals of constant cf av , fi and v. 

Although tangent bifurcations are known to occur generically in chaotic and mixed-dyn- 
amical systems, to our knowledge no such bifurcation has been reported so far in the Henon- 
Heiles system. In |34[ l3*H] we have wrongly surmised that all its periodic orbits existing below 
the barrier energy e = 1 are derivatives of the generic orbits A, B, C (and their repetitions) 
through their bifurcations. This was not correct, as we can demonstrate in the following two 
figures. Here we present a sequence of "spider" -like orbits that are born out of a tangent 
bifurcation occurring at the scaled energy e = 0.988249. Their stability discriminants trM re( j 
are shown in figure EJ and their six genuine shapes in the (y, x) plane in figure El The generic 
pair of a orbits, born with Maslov indices 22 and 23, keeps its shape through three successive 
pitchfork bifurcations at which the orbits c, b and a' are born; these five orbits remain 
hyperbolically unstable at all energies e > 1.016. The latter two bifurcate again, giving birth 
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Figure 4: Stability discriminants trM re( j of the sequence of "spider" orbits in the Henon-Heiles 
potential born out of a tangent bifurcation at the scaled energy e = 0.988249. The numbers give 
their Maslov indices a in the respective stability intervals; the letters refer to the six orbit types 
whose shapes are shown in hgure|51 the horizontal dotted lines correspond to trM re( j = ±2. 

to the orbits c' and a" which remain inverse-hyperbolically unstable at all energies e > 1.016. 
Note that the three orbits a, b, c are reflection-symmetric around the symmetry axes (shown 
by the dotted lines in figure EJ) containing the A orbits, whereas the others are not. The 
Maslov indices, which fulfill again the rules of sectionEU have been obtained with our present 
method. We found, in fact, that the earlier methods of (THj and [TT] could not be applied safely 
here: the use of the intrinsic coordinate systems of these complicated orbits is numerically 
not always stable enough to yield unique results. This actually demonstrates an advantage 
of the present method which works reliably for not too unstable orbits (|trM red | < 40). 
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Figure 5: Shapes of the "spider" sequence of periodic orbits in the Henon-Heiles potential. 
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3.3 A spin-boson system 

We consider a spin-boson system denned by the quantum Hamiltonian 

H = LUb a) a + lo s s 3 + k (a/ + a)(s+ + s_) , s± = si ± is 2 , (45) 

where s a = \a a (a = 1,2,3) are the usual spin operators for s — \ particles {Ti = 1). This 
model has a broad range of applications in atomic, molecular and solid-state physics and in 
quantum optics. In the different fields the Hamiltonian (|45jl bears different names, among 
them "Rabi Hamiltonian" and "molecular polaron model" (see [39J for a review and further 
references) . 

In order to treat this system (semi-) classically, we have to define a phase-space symbol for 
the Hamiltonian (|4"5|). To this purpose we introduce the bosonic operators a = (Qb+iPb)/ V2, 
at = (Qb—iPb) I v2 and take their Wigner transforms to define the canonical bosonic variables 
{Qb, Pb)- For the spin variables n a we use the spin coherent-state symbols of the spin operators 
s a , divided by the value of spin s (see eg jBJ). This leads to the following symbol of the 
Hamiltonian (|43j) 

H = ^ (P 2 + Ql - 1) + ^-n 3 + V2KQ bni . (46) 
The classical equations of motion read 

Q b = u b Pb , Pb = -uj b Qb - \f2nrii , (47) 
n 1 = -uj s n 2 , n 2 = uj s n x - K,Q b n 3 / V2 , n 3 = KQ b n 2 /V2, (48) 

with the constraint n\ + n\ + n\ = 1 . 

We can now introduce the Darboux coordinates (Q s , P s ) by making a stereographic pro- 
jection from the North Pole of the unit n-sphere onto the complex plane and then contracting 
the plane to a disc with radius \/2: 

m = P aA /2-P2-Q2, n 2 = Q ay j2-P?-Q*, n 3 = P 2 S + Q 2 S - 1 . (49) 

Under this mapping, the North Pole is projected onto the boundary of the disc P 2 + Q 2 = 2, 
and the South Pole is projected into the centre of the disc. The Hamiltonian then has the 
form 

H = y (P b 2 + Ql - 1) + y {P 2 S +QI-1) + nQ b P s p{2-P^-Ql). (50) 

This is a two-dimensional harmonic oscillator, perturbed by the nonlinear term proportional 
to k. The representation (j4Uj) is convenient because the equations of motion (|4*7|). (|4*Hj) for 
both boson and spin variables can now be written in a canonical Hamiltonian form: 

Q> = w; p ° = -w: ( ^ M) (51) 

and the equation for the matrizant M(t) is easily found. 

However, as soon as we cross the North Pole on the n-sphere, the equations become 
singular, and one has to switch to the alternative representation 

n 1 = P a y]2-F*-Q*, n 2 = -Q ay /2-P?-Q*, n 3 = 1 - P 2 - Q 2 , (52) 
which corresponds to the projection from the South Pole. 
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The (Qb, Pb) sections of the periodic orbits in this system are similar to those of the 
orbits in the unperturbed harmonic oscillator (« = 0), while the spin components n a on 
the sphere - or, correspondingly, on the (Q s , P s ) disc - evolve substantially with increasing 
coupling constant k. For small k, the two periodic orbits R3 and R5 originate from the South 
and the North Pole, respectively (figured). In the (Qb,Q s ) space they are simple rotations 
with Maslov indices 3 and 5, respectively (figure^, becoming more and more distorted with 
increasing k. At larger values of k they undergo pitchfork bifurcations (see figure OJ), giving 
birth to the orbits P 3 and Q 5 , respectively. The (Qb, Q s ) shapes of the four orbits R 2 , R6, P3 
and Q5 at k = 0.27 are shown in figure |H1 (The subscripts in R CT and Q CT denote again the 
Maslov indices of the respective orbits.) 

We note that in the limit k —>■ the Maslov index of the orbit R3 coincides with the 
Maslov index (a = 3) of the shortest isolated orbit in the unperturbed harmonic oscillator 
(JHHI) with the frequency ratio uj s : Ub = 0.6. The orbit R5 is ill-defined in the representation 
(|4*§|) in this limit, and we have to switch to (p)2"|) instead. For the Hamiltonian (jH3), this 
is equivalent to changing u s — > — to s . Then, the formula (|3^|) yields a = 1 which is equal 
to 5 (mod 4). The difference can be associated with the Maslov index of the matrix which 
transforms the matrizant in the representation ()49|) to the matrizant in the representation 
(}52|) . even though this matrix is not defined at k = 0. 




Figure 6: (Qb,Qs) shapes of orbits R3 (solid 

line) and R5 (dashed line) in the spin-boson Figure 7: Spin components n a of the orbits 

system (f5U|) at ft = 0.05. Other parameters R3 and R5 in the spin-boson system (EH) at 

as in figure |H1 k = 0.05. Other parameters as in figure |H1 
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Figure 8: Shapes Q s versus Qb of the four orbits R2, R6, P3 and Q5 (from left to right) in the 
spin-boson model (|5T?|) at k = 0.27. Other parameter values: ujb = 1.0, u> s = 0.6, E = 4.0. 
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0.0 0.1 0.2 0.3 0.4 

K 

Figure 9: Stability discriminants of the shortest periodic orbits in the spin-boson Hamiltonian (|50j) 
with versus n. Other parameters as in figure El The subscripts denote the Maslov indices. 

The bifurcation scenario is shown in figure El where we plot the discriminants trM re d of 
these orbits versus the parameter k. The orbits R3 and R5 touch the line tr M re d = —2 in the 
stability diagram due to the presence of the discrete reflection symmetry Qb — > —Qb, P s 
—Psi t — > —t. The new orbits P3 and Q5 born at their bifurcations have more complicated 
self-crossing rotational shapes in the (Qb, Q s ) space with a lower discrete symmetry than that 
of R3 and R5 (see figure |BJ). The shape of R3 does not change qualitatively after its successive 
bifurcations, when it becomes R2 and Ri. The same holds for R5 which becomes R6- 

The orbit Q5/Q1 is interesting in the sense that for n = 0.27 it has the Maslov index 
(j = 5, while at k — 0.30 its Maslov index is a = 1. The sign of s is negative in the entire 
interval of existence of this orbit; the change in the Maslov index m is due to a drop of 
the winding number from 3 to 1 near n ~ 0.275. This sudden change of the Maslov index 
by four units, without bifurcation, can be accounted for by a touching of the North Pole 
near k ~ 0.275. It reflects the singularity of the representation ()49)1 and is not felt in the 
semiclassical trace formula (0) where the Maslov enters only modulo multiples of four. 




Figure 10: Spin component s n a of the orbit Figure 11: Spin component s n a of the orbit 
R2 in the spin-boson system (|50| at k = 0.27. R6 in the spin-boson system (|50l) at k = 0.27. 
Other parameters as in figure |H1 Other parameters as in figure |H1 
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The spin components n a of the orbits R2 and Kq ai k = 0.27 are shown in figures ITOl and 
ITTIand those of the orbit Q5/Q1 at k = 0.27 and 0.30, respectively, in figures IT2l and fT3l 




Figure 12: Spin component s n a of the orbit 
Q5 in the spin-boson system (|50l) at k = 0.27. 
Other parameters as in figure |H1 




Figure 13: Spin component s n a of the orbit 
Qi in the spin-boson system ()50l) at k = 0.30. 
Other parameters as in figure |H1 



3.4 Two-dimensional quantum dot with Rashba spin-orbit inter- 
action 

We finally consider a two-dimensional electron gas in a semiconductor heterostructure, lat- 
erally confined to a quantum dot by a harmonic potential. It is modelled by the quantum 
Hamiltonian 

H = -(Pl+f y ) + - (uj 2 x x 2 + uj 2 y f) + 2n{s 2 p x - hPy) ; (53) 

here we put the effective mass of the electrons to be m* = 1. The semi classical treatment 
of this system has been presented recently in [Hj. The classical symbol of the quantum 
Hamiltonian 

H = \{pl+P 2 y ) + \ (u 2 x x 2 + u 2 y 2 ) + k (n 2 p x - n lPy ) (54) 

was considered and the corresponding (semi-)classical equations of motion were studied there. 
Two analytic periodic solutions A x and A y , as well as four numerical solutions D xl , D x2 , 
Dy\, and D^, were found and discussed in jS]. 

We present in table El the Maslov indices of the twelve shortest periodic orbits of this 
system, calculated at the same parameter values as in [Hj. Note that the Hamiltonian (|Hl"J) 
describes a system which is effectively three-dimensional. Therefore, loxodromic blocks occur 
in the monodromy matrix, as well as transitions from a loxodromic block into two elliptic 
blocks without change in the Maslov index ^H] ■ We should also mention that for the calcula- 
tion of the Maslov indices quoted in table El we have used different Darboux representations 
to avoid the problem of crossing the pole of projection. Thus, for D xl and D^ 2 we have 
chosen 

n i = ~Qz\/2-p 2 z -q 2 , n 2 = -{pl + ql ~ 1) , n 3 = p z ^2 - p 2 z - q 2 z , (55) 

which corresponds to the pole of projection located at (0, —1,0), while for D yl and Dy% we 
have projected from the point (1,0,0) 

ni=p 2 z + q 2 z -l, n 2 = q z ^2 - p\ - q 2 z , n 3 = -p z ^2-p 2 z - q 2 z . (56) 
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Hereby we have put q z = Q s and p z = P s . 
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Table 3: St abil ities, Maslov indices and their ingredients of the shortest orbits in the Rashba 
Hamiltonian (|54|) . 



4 Summary 

In this paper we have taken the point of view of practitioners of the semiclassical periodic 
orbit theory. We have formulated a simple calculational recipe for the calculation of Maslov 
indices for isolated periodic orbits that is canonically invariant and does not require the use 
of the orbits' intrinsic coordinate systems. Our work was inspired by two recent formulations 
|13| ITi] which are theoretically very thorough but both have left some practical questions 
unanswered. We have given unique and practicable definitions of the stability angle x an d 
the winding number m, which are the main ingredients of Sugita's formula ((Zj) for the Maslov 
index, and tested them for an integrable and various non-integrable systems. We have found 
that this formula leads to identical results with the method of Wintgen et al [10, and the 
method of Creagh et al JJJ . An alternative definition of stability angle and winding number, 
using a different phase convention, allowed for a direct relation to the decomposition Q given 
in JJJ and lead us to formulate some empirical rules which are useful for the classification 
of periodic orbits in connection with complicated bifurcation scenarios. These rules could 
also be verified in a novel sequence of periodic orbits that we have found in the Henon-Heiles 
system to generate from a tangent bifurcation occurring near the saddle energy. Their shapes 
are so entangled that the use of their intrinsic coordinate systems needed in the methods of 
[TIH E] was numerically not stable enough to yield unique Maslov indices. The present 
method gives unique results as long as these orbits are not too unstable (|trM red | < 40), thus 
demonstrating the practical strength of this method. 

We do not claim to have established any fundamentally new insights here. As a matter 
of fact, some of our steps and observations have been hinted at before in the literature 
jHl EH CHI EE] • Our aim was rather to clarify some practical aspects and to define an easy- 
to-use but canonically invariant method for the calculation of Maslov indices, applicable to 
the most general type of Hamiltonian systems including spin degrees of freedom. We believe 
to have reached this goal and hope that our method turns out to be useful also for other 
practitioners. 
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